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We study the robustness of quantum computers under the influence of errors modelled by strictly 
contractive channels. A channel T is defined to be strictly contractive if, for any pair of density 
operators p,a in its domain, \\Tp — To'\\ 1 < k\\p — for some < k < 1 (here H-^ denotes 
the trace norm). In other words, strictly contractive channels render the states of the computer 
less distinguishable in the sense of quantum detection theory. Starting from the premise that 
all experimental procedures can be carried out with finite precision, we argue that there exists 
a physically meaningful connection between strictly contractive channels and errors in physically 
realizable quantum computers. We show that, in the absence of error correction, sensitivity of 
quantum memories and computers to strictly contractive errors grows exponentially with storage 
time and computation time respectively, and depends only on the constant k and the measurement 
precision. We prove that strict contractivity rules out the possibility of perfect error correction, and 
give an argument that approximate error correction, which covers previous work on fault-tolerant 
quantum computation as a special case, is possible. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Ta. 



I. INTRODUCTION 

Since it was first realized [[l] that maintaining reliable 
operation of a large-scale (multiqubit) quantum com- 
puter in the presence of environmental noise, as well as 
under the combined influence of unavoidable imprecisions 
in state preparation, manipulation, and measurement, 
will pose quite a formidable obstacle to any experimental 
realization of the computer || , many researchers have ex- 
pended a considerable effort devising various schemes for 
"stabilization of quantum information." These schemes 
include, e.g., quantum error-correcting codes (QECC's) 
||, noiseless quantum codes (NQC's) Q, decoherence- 
free subspaces (DFS's) [||, and noiseless subsystems 
(NS's) ||. (The last three of these schemes boil down to 
essentially the same thing, but are arrived at by different 
means.) However, each of these schemes relies for its ef- 
ficacy upon explicit assumptions about the nature of the 
error mechanism. Quantum error-correcting codes || , for 
instance, perform best when different qubits in the com- 
puter are affected by independent errors. On the other 
hand, stabilization strategies that are designed to handle 
collective errors [|]-[| make extensive use of various sym- 
metry arguments in order to demonstrate existence of the 
so-called "noiseless subsystems," i.e., subsystems that are 
effectively decoupled from the environment, even though 
the computer as a whole certainly remains affected by 
errors. 

In a recent publication |Q , Zanardi unified the descrip- 
tion of all above-mentioned schemes via a common al- 
gebraic framework, thereby reducing the conditions for 
efficient stabilization of quantum information to those 
based on symmetry considerations. The validity of this 



framework will ultimately be decided by experiment, but 
it is also quite important to test its applicability in a 
theoretical setting which would make as minimal of an 
assumption as possible concerning the exact nature of 
the error mechanism, and yet would serve as an abstract 
embodiment of the concept of a physically realizable (i.e., 
nonideal) quantum computer. 

In this respect, the assumption of finite precision Q 
of all physically realizable state preparation, manipula- 
tion, and registration procedures is particularly impor- 
tant, and can even be treated as an empirical given. 
This premise is general enough to subsume (a) fundamen- 
tal limitations imposed by the laws of quantum physics 
(e.g., impossibility of reliable discrimination between any 
two density operators with nonorthogonal supports), (b) 
practical constraints imposed by the specific experimen- 
tal setting (e.g., impossibility of synthesizing any quan- 
tum state or any quantum operation with arbitrary pre- 
cision), and (c) environment-induced noise. 

As a rule, imprecisions in preparation and measure- 
ment procedures will give rise to imprecisions in the 
building blocks of the computer (gates) because the preci- 
sion of any experimental characterization of these gates 
will always be affected by the precision of preparation 
and measurement steps involved in such characterization. 
Conversely, precision of quantum gates will affect preci- 
sion of measurements because the closeness of conditional 
probability measures (say, in total variation norm p(|), 
conditioned on the gate used, is bounded above by the 
closeness of any two quantum gates in question Q . 

The central goal of this paper is to offer an argument 
that the concept of a strictly contractive quantum chan- 
nel yields a natural (and very economical) embodiment 
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of the above finite-precision assumption. Defining a suit- 
able distance function d(-, ■) on the set of density opera- 
tors, we say that a channel T is strictly contractive (with 
respect to d) if there exists some k G [0, 1) such that, for 
any pair p, a of density operators, we have the uniform 
estimate 



is addressed in IV D. Finally, in Section |y|, we present 
concluding remarks and outline some open questions and 
directions for future research. 



II. PRELIMINARIES 



d(Tp,Ta) < kd{p,a). 



(1) 



A. On notation 



For instance, the much studied depolarizing channel is 
strictly contractive (with respect to the trace-norm dis- 
tance, to be defined later). If we assume the dominant 
error mechanism of the computer to be strictly contrac- 
tive, then the constant k can be thought of as a quantita- 
tive measure of the computer's (im)precision. While this 
approach may certainly be criticized as reductionist || , 
its merit lies in the fact that it brings out many essen- 
tial features of physically realizable quantum computers 
without invoking more specific assumptions. 

Let us give a "sneak preview" of what is coming up. 
First of all, we establish that the set of all strictly con- 
tractive channels on a particular quantum system (com- 
puter) Q is dense in the set of all channels on Q. Since 
finite-precision measurements cannot distinguish a dense 
subset from its closure jl3| , we draw the conclusion that 
strictly contractive quantum channels (SCQC's) serve 
as a physically meaningful abstract model of errors in 
physically realizable computers. This conclusion is fur- 
ther supported by the fact that, in the presence of 
a strictly contractive error mechanism, the probability 
of correctly discriminating between any two equiprob- 
able quantum states is bounded away from unity (or, 
equivalently, no two density operators in the image of a 
SCQC have orthogonal supports). Next we use a partic- 
ularly important property of SCQC's, namely existence 
and uniqueness of their fixed points, to obtain uniform 
dimension-independent estimates of dccohcrcncc rates of 
noisy quantum memories and computers. We also take 
up the question of possibility of error correction (stabi- 
lization). In this regard, we obtain a rather strong result 
that strictly contractive channels admit no noiseless sub- 
systems. The proof of this claim utilizes ideas from repre- 
sentation theory of operator algebras JT^JT^] and depends 
in an essential way on the property of strict contractivity. 

The paper is organized as follows. In Section || we 
introduce the necessary background on quantum states 
and channels, as well as some relevant facts from oper- 
ator theory. Strictly cont ractive quantum channels are 
introduced in Section 
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where we show that the set 
of strictly contractive channels is dense in the set of all 
channels. Then, in Section IV, we give an interpreta- 
tion of strict contractivity in the framework of optimum 
quantum hypothesis testing (IV A) and then use the fixed 
point theorem for strictly contractive channels to obtain 
estimates on decohe rence rate s of n oisy quantum memo- 
ries and computers (IVB). In IV C , we present the proof 
of nonexistence of noiseless subsystems in the presence of 
SCQC's. The possibility of approximate error correction 



In this paper we will adhere to the following notational 
conventions. First of all, the operator adjoint to X will 
be denoted by X*, as is usually done in mathematical 
physics literature. Secondly, density operators will be de- 
noted by p and a (with subscripts, whenever necessary). 
We will use capital Latin letters to denote all other op- 
erators; whenever no ambiguity may arise, the action of 
a mapping X on a density operator p will be written as 
Xp. Finally, the Pauli matrices will be written as q, 
% G {1,2,3}. 



B. States 

Let Ti. be a finite-dimensional Hilbert space associated 
with the computer Q, and let B(Ti) be the algebra of all 
bounded operators on TL [since dimH < oo, the quali- 
fication "bounded" is patently unnecessary, but we will 
retain the notation B(H), following standard usage]. The 
set 



S(H) := {pe B(H)\p>Q;trp= 1} 



(2) 



is the set of all density operators (states) of Q. We can 
define a few norms on B(Ti.); since H. is finite-dimensional, 
all norm topologies on it arc equivalent. First, we have 
the operator norm 



\X\ 



sup 



VX g B(H). 



(3) 



We can also define the class of Schatten p-norms |T6[ . For 
any X G B{H), we let \X\ := (X*X) 1 / 2 , so that 



\X\\ p := (tr \X\ P ) 1/P , 



\/X g B(H);p = 1,2, 



(4) 



The Schatten 1-norm is better known as the trace norm; 
in the case p — 2, we recover the Hilbert-Schmidt norm. 
In fact, for any X G B(H), \\X\\ p -> ||X|| as p -> oo. For 
this reason, we can identify the operator norm ||-|| with 
11*11 oo ■ ^ u th csc norms possess a very important property 
of unitary invariance pi : for any unitaries U, V and any 
X G B(H), we have 



\uxv\\ p = 



\x\ 



P=l,2, 



(5) 



The trace norm can be given a natural interpretation as 
a distance between density operators . First of all, for 
any p G S(Ti), we have \\p\\i = 1- Of especial importance 
is the fact that, for any pair p, a G S(H), the trace-norm 
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distance \\p — o~\\-, achieves its maximum value of 2 if and 
only if pa = (i.e., if and only if p and a have orthogo- 
nal supports). In the case of two pure states \c/>){^>\ and 
|^>)(0|, this condition reduces to (<f>\ip) = 0, i.e., the cor- 
responding state vectors <j>, ip G H are orthogonal. As we 
will see in Section IV, the trace- norm distance also fig- 
ures prominently in the framework of optimal quantum 
hypothesis testing. 

To close our discussion of states, we give two impor- 
tant characterizations of the trace-norm distance. Let 
X be a selfadjoint operator. Then we can write X as a 
difference of two positive operators with orthogonal sup- 
ports: X = X+-X-, where X± := (\X\±X)/2. This is 
referred to as the orthogonal decomposition of X. Then 



\X = X, 



\x\ 



X- , and 



IX, 



\x\ 



tr X+ + tr A_ 



Now let p, a be a pair of density operators. 
p — a = R + — i?_ , we observe that tr (p — a) 
plies tr R + = tr i?_ , and hence 

||p - = 2tri?+ . 

Another useful relation is 

lip — er||, = 2 max tr F(p — er), 
1 0<F<I 



(6) 

Writing 
= im- 



(7) 



(8) 



where the inequality < F < 1 should be taken to mean 
F > and 1 — F > 0. In fact, since the set of all such F is 
convex, the maximum of the linear functional in Eq. (0) is 
attained on an extreme point of {F | < F < 1}, namely 
on the projector P + defined by 



P+R- = 0. 



(9) 



C. Channels 

In quantum theory, the reversible evolutions of a 
closed quantum system correspond to the automorphisms 
S(H) -> US{H)U* with a unitary U, i.e., for any 
p G 3(H), we have p i— > Tup := UpU*. The map 
T[j : S(H) — > S(H) is affine, trace-preserving, and posi- 
tive (we will call a map positive if it takes positive opera- 
tors to positive operators). Since any affine map on den- 
sity operators can be uniquely extended to a linear map 
on selfadjoint operators |Q, we can take T\j to be linear. 
Most importantly, Tjj is invertible with T u ~ 1 p := U*pU. 

The general irreversible evolution T : S(TL) — > S(TL) 
of an open quantum system will no longer be given by 
an automorphism U ■ U*. We must accordingly modify 
the requirements imposed on T. It is obvious that we 
have to drop the invertibility condition, so that T is now 
a trace-preserving positive linear map on S(Tl) (we can 
then extend it, by linearity, to selfadjoint trace-class op- 
erators). However, positivity alone is not sufficient. In 
order for T to represent a physically admissible evolution, 



it must be completely positive |12|, i.e., the map T <S> id n , 
where id„ is the identity operator on the space M„(C) 
of n x n complex matrices, must be positive for all n. 
Unitary evolutions Tjj obviously satisfy all these require- 
ments. In fact, we can also include such transformations 
as measurements into this framework by requiring all ad- 
missible evolutions to be trace-nonincreasing completely 
positive linear maps on S(Tl), i.e., for any p G S(H), we 
have tr T p < tr p, so that 



tvTp' 



(10) 



Then tr Tp can be naturally interpreted as the condi- 
tional probability of transformation T occurring given 
that the system is initially in the state p. We will call 
any trace-preserving completely positive linear map on 
S(TL) a channel. 

There are many useful structure theorems for com- 
pletely positive ma ps. For instance, the Kraus repre- 
sentation theorem |L2| states that, for any completely 
positive map T : B(H) — ► B(H), there exists a collection 
{Ki] of bounded operators such that 



T(X) = K l XK*, VA G BiH). 



(11) 



If T is trace-nonincreasing, then we have the bound 



E 



KtKi < I, 



(12) 



where equality is achieved if and only if T is trace- 
preserving. In addition, if T is a unital channel, i.e., 
T(l) = 1, then we also have 



J2K t K* = l. 



(13) 



Now we must adopt a suitable metric on the set of all 
completely positive maps on B(H,). One possible candi- 
date is the metric induced by the operator norm, 



|T|| := sup UrpOU 

X6B(W);||X||=1 



(14) 



Unfortunately, the operator norm is rather ill-behaved 
p9| : it is not stable with respect to tensor products. In 
particular, if T is a completely positive map, then the 
norm ||T(g>id n || can in general increase with n. A good 
choice then is the metric induced by the norm of complete 
boundedness Jl9| (or cb-norm), defined as 



\T\ 



cb 



sup ||T ® id„ 



(15) 



This norm has appeared, under different guises, in 
Refs. p7| ], ]20[ |, and |21|. For any selfadjoint trace-class 
operator X on 7i and any two maps S, T on B{TL) with 
finite cb-norm (in the case of finite-dimensional 7i, this 
is always true [pi), we have the relations [pol 
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ITOIIx <\\T\\ ch ll^llx, 
||T5|| cb <||T|| cb ||5|| cbJ 
|r®S|| cb =||T|| cb ||5|| cb . 



(16) 
(17) 
(18) 



|TH eb = l. 



Furthermore, for any channel T, we have [E2| 

If two channels T, S are close in cb-norm, then, for any 
density operator p, the corresponding states Tp, Sp are 
close in trace norm since, from Eq. (|l6|), it follows that 



\\Tp-Sp\\ 1 = \\(T-S)p\\ 1 <\\T-S\ 



cb • 



(19) 



In fact, the above estimate cannot be loosened by adjoin- 
ing a second system with Hilbert space K, in some state 
a, entangling the two systems through some channel K 
on S(H (g> /C), and then comparing the channels T (g> R 
and S (£> R, where R is some suitably chosen channel on 
S(IC). This is evident from the estimate 

|| (T ® R)K (p ® a) -(S®R)K(p Qa)^ < ||T-S|| cb , 

(20) 

which can be easily obtained by repeated application of 
Eqs. (|l6|)-(|is|). In other words, as far as the cb-norm 
distinguishability criterion is concerned, entangling the 
system with an auxiliary system will not improve dis- 
tinguishability of the channels T and S. The cb-norm, 
however, is an extremely strong distinguishability mea- 
sure: its definition already accounts for optimization with 
respect to entanglement and input states over Hilbert 
spaces of very large (but finite) dimension. There exist 
weaker measures of channel distinguishability (such as 
the channel fidelity p3fl), which describe how channels 
may be distinguished with only finite resources. Using 
these weaker criteria, one may show that entanglement 
does improve practical distinguishability of both states 
and channels H). 

Before we go on, we must mention that, for the present 
purposes, we only need to consider channels that map op- 
erators on some Hilbert space H to operators on the same 
Hilbert space. In general, this does not have to be true. 
For instance, if the Hilbert space in question is a tensor 
product Tli ® TC2, then the partial trace over TL2 can be 
treated as a channel ti2 ■ S(Hi <8> H2) — » S(Hi). 



D. Some facts from operator theory 

We close Section || by listing some facts from op- 
erator theory, which will be necessary in the sequel. 
Let X be a metric space with the corresponding met- 
ric d(-, •). An operator A : X — > X is called a contrac- 
tion if, for any x,y £ X, d{Ax,Ay) < d(x,y), and a 
strict contraction if there exists some k £ [0, 1) such that 
d(Ax,Ay) < kd(x,y). If X is a complete metric space, 
then the contraction mapping principle |^5| states that 
any strict contraction A on X has a unique fixed point. 
In other words, the problem Ax = x has a unique solu- 
tion on X. If y is a closed subset of X, then it follows 



that any strict contraction A : y — > y has a unique fixed 
point on y. 

Strict contractivity is a remarkably strong property. 
Indeed, if we pick any y £ y, then the sequence of iterates 
A n y converges to the fixed point yo of A exponentially 
fast, because 



d(A n y,y ) = d(A n y,A n y ) < k n d(y,y ). 



(21) 



This fact is of tremendous use in numerical analysis when 
one wants to solve a fixed-point problem Ay — y via it- 
eration method with some initial guess y. If the operator 
A is a strict contraction on a closed subset of a com- 
plete metric space, then, for any choice of y, the itera- 
tion method is guaranteed to zero in on the solution in 
OQoge^ 1 ) steps, where e is the desired precision. 

It should be noted that existence and uniqueness of a 
fixed point of some operator A are, by themselves, not 
sufficient to guarantee convergence of the sequence of it- 
erates A n y for any point y in the domain of A. Indeed, 
according to the Leray-Schauder-Tychonoff theorem p5[ , 
any continuous map on a compact convex subset of a lo- 
cally convex space X has at least one fixed point. Fur- 
thermore, any weak contraction on a compact subset C of 
a Banach space, i.e., a map W : C — > C with the property 
\\Wx — Wy\\ < \\x — y\\ for any x,y 6 C, has a unique 
fixed point (^6|. The key to the rapid convergence in 
Eq. (|2l]) is the fact that a strict contraction A : y — > y 
shrinks distances between points of y uniformly. 



III. STRICTLY CONTRACTIVE QUANTUM 
CHANNELS 

A. Definiton and examples 

Let H be the finite-dimensional Hilbert space associ- 
ated with some quantum system Q. Then, as follows eas- 
ily from Eq. (ttq), any channel T on S(Ti) is a contraction: 



(22) 



In other words, no channel can make any p, a more dis- 
tinguishable. For a channel T, we define the contractivity 
modulus 



k(T) := sup 



\\Tp-Ta\l, 



P,<7£S(H) WP-O'Wl 



(23) 



Any channel T with k(T) < 1 is strictly contractive, and 
thus has a unique fixed point px 6 S(TL). 

The depolarizing channel D p , < p < 1, whose action 
on an arbitrary p £ S(Ti) is given by 



D p p :=P^ + (1 -p)p, 



(24) 



where d — dim7i, is manifestly strictly contractive with 
k(D p ) = 1 — p. The maximally mixed state 1/d is the 
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unique fixed point of D p for any p. When d = 2, so that 
H = C 2 , the action of D p on S(H) can be visualized as a 
uniform rescaling of the Bloch-Poincare ball by a factor 
of k(D p ), and the term "strictly contractive" thus be- 
comes especially apt. It also turns out that, for any two 
depolarizing channels D p and D q , their tensor product is 
also strictly contractive. In order to show this, we use the 
fact that a density operator on C 2 (g) C 2 can be written 



P = 7 




ft 



(25) 



A./ 



where the vectors a := (ai, a 2 , 0:3) and (3 := /3 2 , P3) 
are the coherence vectors of the first and second qubit 
respectively, while the matrix 8 with entries Oki is called 
the correlation tensor of p. The action of the depolariz- 
ing channel on an arbitrary operator X G B(C 2 ) can be 
described as 



D p (X) 



p(trX)~ 



(l-p)X 



[if X is a density operator, this reduces to Eq. (p4j)]. 
•Dp ft = (1 - p)ft, which yields 

(D p <g> = i f 1 <g) 1 + (1 -p) ^ a fc 9 fe ® 1 
+ (l-p)(l-«)5^fl w ?*<8)« 



(26) 



Thus 



(27) 



A./ 



It is then straightforward to verify that T„, 



Dr,®D n 



is strictly contractive with K,(T pq ) — max[(l — p), (1 — q)]. 
In particular, the channel D p (g _D p is strictly contractive 
with k(-D p (g Dp) = k(D p ) = I— p. Strict contractivity of 
the product channel D p (g D p provides an alternate ex- 
planation of the fact that the use of entanglement cannot 
improve distinguishability of classical signals transmitted 
through the depolarizing channel p8fl . 

In fact, since the trace-norm distance between any two 
density operators on C 2 is just the Euclidean distance 
between their Bloch-Poincare vectors, any strictly con- 
tractive channel on 5(C 2 ) can be pictured as a rescaling 
of the Bloch-Poincare ball (which may not be isotropic, 
as long as the maximum of the scaling ratio over all di- 
rections is strictly less than one), possibly followed by 
translation and rotation. As shown in Ref. p9| , for any 
channel T on £>(C 2 ) [which is just the space M2(C) of 
2x2 complex matrices] , there exist unitaries U, V and 
vectors v.tel 3 such that 



Tp = U[T v , t {VpV*)]U\ 



(28) 



where the action of T v ,t is defined, with respect to the 
basis {1, 91, 92, 93}, as 

T V: t(^ol + w-9) :=w !+[t + (diagv)w] .9. (29) 

Assuming that v and t are such that the map T is in- 
deed a channel p9]|, we see that T is strictly contractive 
whenever max^^ 2 3 j < 1. In fact, the contractivity 
modulus of T satisfies 



k(T) — max \vj 
»e{i,2,3}' 



(30) 



If T\ , T 2 are unital strictly contractive channels on 
S(C 2 ), then the product channel T\ ® T 2 is strictly con- 
tractive on S(C 2 (g C 2 ). Given a representation ( |2q ) of a 
channel T on 6>(C 2 ), we see that T is unital if and only 
if t = 0. Specifically, for T\ and T 2 we have 



T 1 p=U 1 {T Vu0 (V lP V 1 *)}U?, 
T 2 p^U 2 [T V2fi (V 2P V 2 *)}U;. 



(31) 
(32) 



Thus the action of the channel T\ <g T 2 on a density 
operator p over C 2 ® C 2 is a successive application of 
the unitary channel (V\ ®V 2 ) ■ (Vf <E)V 2 ), the rescal- 
ing transformation T V1 .0 <E> T V2 -0 , and the unitary chan- 
nel (U\ (g U 2 ) ■ (U* (g U 2 ) to p. By unitary invariance 
of the trace norm, we only need to consider the effect 
of Tvi.o ® T V2 ^ . Writing p in the form of Eq. (p5|), we 
obtain 



{T Vu0 <g> T V2 . )p ■■ 



ft 



y^fc flfcfft <g 9; 



(33) 



A',/ 



where uj ,i S {1,2}, j G {1,2,3}, denotes the jth com- 
ponent of Vj. By inspection, 



k(Ti <g T 2 ) 



max 



,W 



(34) 



If at least one of the channels Ti and T 2 is not unital, the 
tensor product channel T — 7\ (g T 2 may not be strictly 
contractive, even if Ti and T 2 are. This stems from the 
fact that, in this case, the effect of T on the correlation 
tensor O of p is determined not only by T, but also by p 
through the coherence vectors a and (3. 

It is quite easy to see that any unital strictly con- 
tractive channel maps all states to mixed states. Let 
d = dim7Y. Then, for any unit vector ip G Ti., we have 



MM 



2(d- 1) 



(35) 
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(this can be readily proved by expanding 1 with respect 
to an orthonormal basis containing ip). Now suppose 
that T is a strictly contractive unital channel that maps 
\4>){ip\ to some other pure state \<fr)(<f>\- Then 



Tl 

TWW- — 



\m\-- d 



Furthermore, we must also have 



Tl 



<k(T) 



2(d-l) 



(36) 



(37) 



Hence, k(T) > 1, which is a contradiction, since k(T) < 1 
for any strictly contractive channel. 

We can show that there exist channels that are not 
strictly contractive, and yet contain no pure states in 
their image. Let T be an arbitrary channel on S(TL). 
According to the Leray-Schauder-Tychonoff theorem, T 
has at least one fixed point on S(H). Let us adjoin an- 
other system with the associated Hilbert space K,. Then 
the channel T £g> id on S(TL <g> JC) cannot be strictly con- 
tractive because, for any fixed point px of T and any 
a G S(1C), the product density operator p? <8> o is a fixed 
point of T £g> id. The channel T <£> id is not even weakly 
contractive, because it preserves the trace-norm distance 
between any two of its fixed points. However, if the im- 
age of <S(H) under T contains no pure states, then the 
image of S(H ® K.) under T <E> id contains no pure states 
either because of the relation J3(| 

P Mk) S{Tp)= P esf^^^' (38) 
where S(p) is the von Neumann entropy of the state p. 



B. Strictly contractive channels are dense in the set 
of all channels 

As we have mentioned in the Introduction, finite- 
precision measurements cannot distinguish a dense sub- 
set from its closure. Let us make this statement more 
precise. Suppose we are presented with some quantum 
system Q in an unknown state p, and we are trying to 
estimate the state. Any physically realizable apparatus 
will have finite resolution e, so that all states p' with 
1 1 yO — f/\\, < e are considered indistinguishable from p. 
Now, if TL is the Hilbert space associated with Q, and 
if £ is a dense subset of S(H), then, by definition of a 
dense subset, for any e > and any p G S(7i), there will 
always be some a G S such that \\p — cr||j < e. 

The same reasoning also applies to distinguishability 
of quantum channels, except now the appropriate mea- 
sure of closeness is furnished by the cb-norm. Thus, if 
an experiment utilizes some apparatus with resolution e, 
then any two channels T, S with ||T — S\\ ch < e are con- 
sidered indistinguishable from each other. There is, how- 
ever, no fundamental difference between distinguishabil- 
ity of states and channels because any experiment pur- 
porting to distinguish any two channels T and S consists 



in preparing the apparatus in some state p and then mak- 
ing some measurements that would tell the states Tp and 
Sp apart from each other. Then, since for any state p, 
\\Tp — Sp\\ l < \\T — S\\ ch , the resolving power of the ap- 
paratus that will distinguish between T and S is limited 
by the resolving power of the apparatus that will distin- 
guish between Tp and Sp. 

In this regard, we have the following 
Proposition 1. Let C(H) be the set of all channels on 
S(H), where TL is the Hilbert space associated with the 
system Q. Then the set C SC (H) of all strictly contractive 
channels on S(TL) is a ||-|j cb -dense convex subset of C(TL). 
Proof. We show convexity first. Suppose T\,T 2 G 
C SC {H). Form the channel S := ATi + (1 - A)T 2 , 
< A < 1. Then, for any p, a £ S(H), we have the 
estimate 

\\S P - SaW, < A \\T lP - TM\i + (1 - A) \\T 2 p - T 2 a\\ x 
< [A K (T 1 ) + (l-A)«(T 2 )]||p-a|| 1 . (39) 

Defining k := max [/c(Ti), k(T2)], we get 



WSp-SaWi < K\\p-a\\ l 



(40) 



Since T\,T 2 are strictly contractive, k < 1, and therefore 
S G C SC (H). To prove density, let us fix some a € S(H). 
Now the map K a : p G S{TL) i — >■ cr is obviously a channel, 
which is furthermore trivially strictly contractive because 
it maps all density operators p to a. Given e > 0, pick 
some positive n such that 1/n < e. For any T G C(H), 
define 



T n := — K„ 
2n 



1 - — ] T. 

2n 



(41) 



Clearly, T n G C SC (H), and the estimate 

r-^|| cb =^]]T-^|| cb <i<e (42) 

finishes the proof. ■ 
This proposition indicates that, as far as physically re- 
alizable (finite-precision) measurements go, there is no 
way to distinguish any channel T from some strictly con- 
tractive T' with || T — T"|| cb < e, where e is the resolution 
of the measuring apparatus. In this regard, it is interest- 
ing to mention that any channel T with ||T — id|| b < e 
(for some sufficiently small e > 0) cannot be distin- 
guished from a depolarizing channel. Indeed, let M be 
the channel that maps all density operators p to the maxi- 
mally mixed state where d = dimTL. Then it suffices 
to pick some 



n > 



|M-id| 



-r-idn cb 



(43) 



so that 



G 



\\T - D 1/n \\ ch < \\T - id|| cb + (1/n) ||M - id|| cb < e. 

(44) 

We note that a convex combination of any channel 
with a strictly contractive channel is a strictly contractive 
channel. Let T G C be an arbitrary channel, and sup- 
pose that T' 6 C sc [from now on, we will not mention the 
Hilbert space 7i when talking about channels on S(H), 
unless this omission might cause ambiguity] . Define, for 
some < A < 1, the channel S := AT + (1- A)T'. Then 



\\S P - Sa^ < A \\Tp - Tal^ + (1 - A) \\T'p - T'a 
< [X + il-XMT'^Wp-aW,. 



(45) 



Since A + (1 - A)k(T') < 1, we conclude that S G C sc . 

Finally, we mention that Proposition 1 implies that 
the set C* c of all unital strictly contractive channels is a 
dense convex subset of the set C 1 of all unital channels. 



IV. IMPLICATIONS FOR QUANTUM 
INFORMATION PROCESSING 



A. Optimum quantum decision strategies 

In this subsection we explore an interesting connec- 
tion between the contractivity modulus of a channel and 
quantum detection theory |3l| . The archetypal problem 
in quantum detection theory is that of optimum M-ary 
detection. A quantum system is prepared in a state p, 
drawn from a collection {p^f^ of M density operators, 
where pi is selected with probability 7Tj. Our task is to 
determine, as accurately as possible, which state pi has 
been drawn. On this system we can perform a measure- 
ment described by a positive operator-valued measure 
(POVM), i.e., a collection {Ti}^!, of operators that sat- 
isfy 



< Fi < 1, 

M 

5> = 1. 



1=1, 



,M, 



(46) 
(47) 



We seek a POVM that would solve the optimization prob- 
lem 



P c = max 2J TTi tr FiPi 



(48) 



subject to the constraints ( flq ) and (f47|). The quantity 
being maximized in j4q ) is the probability of correct de- 
cision using the POVM {Fi}. We will only consider the 
case M = 2, wherein the system can be in the state p\ 
with probability m, or in the state pi with probability 
7T2 = 1 — 7Ti. In this case, we are considering two-element 
POVM's {F, 1 - F} with < F < I, and the optimiza- 
tion problem (E|) takes the form 



max \wi tr F p\ +TT2 tr (1 
0<F <i ' 



F) P 2 



or, equivalently, 



(49) 



(50) 



We interpret trFpi as the conditional probability that 
the measurement using the POVM {F, 1 — F} correctly 
determines the state of the system to be p\\ similarly, 
tr (1 — F)p2 is the conditional probability that the state 
P2 is identified correctly. Then X,vFp2 and tr(l — F)p\ 
respectively are the conditional probabilities of mistaking 
P2 for p\ and vice versa. 
We can easily show that 



1 1 „ 



K2Pl\ 



(51) 



Writing down the orthogonal decomposition ~K\P\ — 
TT2P2 = R+ — R- ) we get tr R+ = n\ — n2 + tr i?_ . Now 



max trF(i? + - R) = tri?_, 

0<F<1 



(52) 



where the maximum is attained by choosing F to be 
a projection operator with FR + = R + and = 0. 

Since 



tr \ttiPi — t^2Pi\ = tri?+ + tri?_ = 2 tr 



-7T 2 



7Tl, 

(53) 



we finally arrive at Eq. (pl|), which clearly exhibits the 
role of the trace-norm distance in optimum quantum hy- 
pothesis testing. It can be proved B2| that P c = 1 if 
and only if p\p2 — 0, in which case HT 7r i/ 5 i — 7r 2P2|li = 

TTl + 7T 2 = 1. 

Now suppose that the state of the system is given by 
one of two equiprobable density operators p\,p2- Sup- 
pose furthermore that p\P2 — 0, so that \\p± — paid = 2. 
Then there exists a measurement that would correctly 
distinguish between p\ and p2 with probability one. Since 
any channel T will generally decrease the trace-norm dis- 
tance ||/9i — P2II1, it can happen that the states Tpi and 
Tp2 no longer have orthogonal supports, and thus the op- 
timum decision strategy will fail with nonzero probability 

Pe = 1 - Pc- 

If T is a weakly contractive channel, then no two den- 
sity operators in its image have orthogonal supports, but 
the probability of error P e can, in principle, be made ar- 
bitrarily small. If, however, T is strictly contractive, then 
we have the trivial, but important, 

Lemma 1. Let T be a strictly contractive channel. 
Then, for any pair p± , P2 of equiprobable density oper- 
ators, the optimum decision strategy for Tp\ and Tp2 is 
such that 



Pr < 



1 + k(T) 



< 1. 



(54) 
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Proof. Fix a pair pi,p2 of density operators. Then, 
using Eq. (|5l|), we get 



11., ,1 k(T) 



\\P1-P2\\i- (55) 



Since \\pi - p 2 \\ 1 < 2, we obtain Eq. (pi]). ■ 
We note that the statement of the above lemma can 
be extended to general channels. For instance, if T is a 
channel with the property that there exists at least one 
pair p, a of density operators such that \\Tp~ Ta\\ 1 — 
\\p — a\\ lt then the bound ( |5~i| ) is obviously 



Pr < I- 



(56) 



If T is a weakly contractive channel, then the inequal- 
ity ( p6| ) becomes strict, but P c can, at least in principle, 
be made arbitrarily close to one. This is decidedly not 
the case for a strictly contractive channel, in which case 
Lemma 1 states that for any pair of equiprobable density 
operators p, a, the probability P c of correctly discrimi- 
nating between them is bounded away from one. 

The discussion in this subsection lends further support 
to our argument that strictly contractive channels serve 
as an abstract model of errors in physically realizable 
quantum computers. In any realistic setting, no event oc- 
curs with probability exactly equal to unity. For instance, 
we can never prepare a pure state \ip)(ip\, but rather a 
mixture (1 — e)\ip){ip\ + ep, where both e and p depend on 
the particulars of the preparation procedure. Similarly, 
the measuring device that would ideally identify IV'XV'I 
perfectly will instead be realized by (1 — 5)\?p)(ip\ + SF, 
where S and the operator F, < F < 1, are again deter- 
mined by practice. If we assume that, in any physically 
realizable computer, all state preparation, manipulation, 
and registration procedures can be carried out with fi- 
nite precision, then it is reasonable to expect that there 
exist strict bounds on all probabilities that figure in the 
description of the computer's operation. 



B. Decoherence rates of noisy quantum memories 
and computers 

So far, we have established two important properties of 
strictly contractive channels. Firstly, any channel T can 
be approximated, in cb-norm, by a strictly contractive 
channel T", and there will always be some finite-precision 
measurement which will not be able to distinguish T from 
T' . Secondly, any measurement that would, in principle, 
distinguish some pair /?, a of density operators with cer- 
tainty, will fail with probability at least [1 — k(T)}/2 in 
the presence of a strictly contractive error channel T. 
The latter statement can also be phrased as follows: no 
two density operators in the image TS(7i) of S(H) under 
some T 6 C sc have orthogonal supports; furthermore, the 
trace-norm distance between any two density operators 
in TS{TL) is bounded from above by 2k(T). 



In this subsection, we obtain dimension-independent 
estimates on decoherence rates of quantum memories 
and computers under the influence of strictly contractive 
noise and without any error correcti on (th e po ssibili ty of 
error correction will be addressed in IV C and IV E ) . 

We treat quantum memories (registers) first. Suppose 
that we want to store some state po S S(Tl) for time t in 
the presence of errors modelled by some strictly contrac- 
tive channel T . Let r be the decoherence timescale, with 
t Ct, and let n = \t/r~\. The final state of the register 
is then p n — T n P Q. If P t is the unique fixed point of T, 
then 

\\ Pn - prW, = ||T> - T n p T \\ x < K{T) n ||po - p T \\ x . 

(57) 

In other words, the state po, stored in a quantum regis- 
ter in the presence of strictly contractive noise T, evolves 
to the fixed state P t of T, and the convergence is in- 
credibly rapid. Let us consider a numerical example. 
Suppose that k(T) — 0.9, and that initially the states 
po and pr have orthogonal supports, so ||po — Pt\\\ = 2. 
Then, after n = 10 iterations (i.e., t = lOr), we have 
\\Pn — Pt\\i < 0.697, and the probability of correct dis- 
crimination between p n and px is only 0.674. Note that 
the decoherence rate estimate 



r(n;p,T) 



\\Pn - grlli 
\\Po -Prili 



< K(T) n 



(58) 



does not depend on the dimension of 7i, but only on the 
contractivity modulus k(T) and on the relative storage 
duration n. In other words, quantum registers of any size 
are equally sensitive to strictly contractive errors. 

Obtaining estimates on decoherence rates of comput- 
ers is not so simple because, in general, the sequence 
{p n }, where p n is the overall state of the computer after 
n computational steps, does not have to be convergent. 
Let us first fix the model of a quantum computer. We 
define fl33| an ideal quantum circuit of size n to be an 
ordered n-tuple of unitaries U i: where each Ui is a ten- 
sor product of elements of some set Q of universal gates 
|?4|| . The set Q will in, general, be a dense subgroup 
of the group U(H) of all unitary operators on TL. For 
some error channel T, a T -noisy quantum circuit of size 
n with k error locations is an ordered (n + fc)-tuple con- 
taining n channels Ui :— Ui ■ U* , where the unitaries Ui 
are of the form described above, as well as k instances 
of T. We will assume, for simplicity, that each T is pre- 
ceded and followed by some Ui and Ui+\ respectively. 
Based on this definition, the "noisiest" computer for fixed 
T and n is modelled by a T-noisy quantum circuit of 
size n with n error locations, i.e., a 2n-tuple of the form 
(t/i, T, U 2 , T, . . . , U n , T). If the initial state of the com- 
puter is po, then we will use the notation 



WTUApo 



(59) 
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to signify the state of the computer after n compu- 
tational steps. In the above expression, the product 
sign should be understood in the sense of composition 
Totl n o...oToU 1 . 

Given an arbitrary sequence of computational steps, 
the sequence {p n }, defined by Eq. (|59|) (assuming that n 
is suficiently large, i.e., the computation is sufficiently 
long) need not be convergent. However, if the chan- 
nel T is strictly contractive, then, for any e > 0, there 
exists some Nq such that, for any pair of initial states 
Po, <J G S(TL), the states p n , o~ n , n > N 0l will be indistin- 
guishable from each other. In other words, any two suf- 
ficiently lengthy computations will yield nearly the same 
final state. Using Eq. (|59|), as well as unitary invariance 
of the trace norm, we obtain 



On-(?n\\i= \Y\ T Uij (po - oo) 

</cCr) n llpo-0-o|li. 



(60) 



Now suppose that, at the end of the computation, we 
perform a measurement with precision e, i.e., any two 
states p,a with \\p — < e are considered indistin- 
guishable. Then, if the computation takes at least No = 
[log(e/2)/ logK(T)] steps, we will have ||p n — c n |li < e 
for all n > Nq. For a numerical illustration, we take 
k(T) = 0.9 and e = 0.01, which yields N = 50. In 
other words, the result of any computation that takes 
more than 50 steps in the presence of a strictly contrac- 
tive channel T with k(T) = 0.9 is untrustworthy since we 
will not be able to distinguish between any two states p 
and a with ||p — (t||j < 0.01. Again, No depends only on 
the contractivity modulus of T and on the measurement 
precision e, not on the dimension of TL, at least not explic- 
itly. We note that, if the state of the computer is a den- 
sity operator over a 2 fc -dimensional Hilbert space, then 
any efficient quantum computation will take 0(Poly(/c)) 
steps, and therefore the sensitivity of the computer to 
errors grows exponentially with k. 

Let us consider some cases where the sequence {p n } 
does converge. Suppose first that the channel T E C sc 
is unital. Then, since each channel Ui is unital as well, 
the sequence {p n } converges exponentially fast to the 
maximally mixed state l/d, where d = diraTL. If the 
computation employs a static algorithm, i.e., Ui — lJ for 
all i (this is true, e.g., in the case of Grover's search al- 
gorithm |35|), then the channel S :— TV is also strictly 
contractive, and k(S) = k{T) by unitary invariance of 
the trace norm. Denoting the fixed point of S by ps, we 
then have 

\\Pn -psh= \\S n P0 - SysW, < n(T) n \\ Pn -psW,, 

(61) 

i.e., the output state of any sufficiently lengthy compu- 
tation with a static algorithm will be indistinguishable 
from the fixed point ps of S = TU. 



C. Impossibility of perfect error correction 

After we have seen in the previous subsection that 
quantum memories and computers are ultrasensitive to 
errors modelled by strictly contractive channels, we must 
address the issue of error correction (stabilization of 
quantum information). Since we have not made any 
specific assumptions (beyond strict contractivity) about 
the errors affecting the computer, it is especially impor- 
tant to investigate the possibility of error correction, if 
only to determine the limitations on the robustness of 
physically realizable quantum computers from the foun- 
dational standpoint. 

First of all, strict contractivity rules out the possibil- 
ity of perfect quantum error-correcting codes Let 
us recall the basics of QECC's. We seek to protect a 
quantum system with a /c-dimensional Hilbert TL space 
by realizing it as a subspace JC (called the code) of a 
larger n-dimensinal Hilbert space TL Cl known as the cod- 
ing space. In other words, the Hilbert space TL is embed- 
ded in the coding space TL C via the isometric encoding 
operator V^nc : TL — > /C. Now, for any channel T on 
S(TC C ), a theorem of Knill and Laflamme H asserts that 
a subspace /C of TL C can serve as a T-corrrecting code if 
and only if there exists some channel S on S(TL C ) such 
that ST\ K = id. In other words, S is the left inverse of 
the restriction of T to 5(/C). However, if the channel T 
on S(TL C ) is strictly contractive, then no subspace K, of 
TL C is a T-correcting code. Suppose, to the contrary, that 
such a subspace /C exists, and let {e^} be any orthonor- 
mal basis of JC. Then there also exists some channel S 
on S(TL C ) that satisfies the Knill-Laflamme condition for 
T and /C. Thus 

\\ST(\e,)(e,\ - \e v ){e v \)\\i = IIMM " l^^lllj (62) 

for all p. v. But, using Eq. ( |l6|) and strict contractivity 
of T, we also have 



MM ~ |e I ,)(e„||| 1 < k(T) |||e M )(e A1 | - |e„)(e„||| 1 



(63) 



which yields k(T) = 1. Since T is strictly contractive, 
this is a contradiction, and therefore no subspace fC of 
TL C is a T-correcting code. 

Before we go on, we must mention that the Knill- 
Laflamme theorem provides also for approximately cor- 
rectable channels. That is, let {K{\ be the set of the 
Kraus operators of some channel T on S(TL C ). For any 
subset A of {Ki}, we can define the completely positive 
map T\ via 



T A {X) := KiXKt, VA e B(TL C ). 



(64) 



KiEA 



Then a subpace IC of TL C can serve as a TA-correcting 
code if there exists some channel S on S(TL C ) such that 
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ST A \ K oc id. If \\T-T A \\ ch is sufficiently small, then 
the errors modelled by the channel T are approximately 
correctable. Thus, in and of itself, the impossibility of 
perfect error correction for strictly contractive channels 
is not likely to be a serious problem. 

However, strict contractivity also proscribes the ex- 
istence of noiseless subsystems in the sense of Knill- 
Laflamme- Viola || and Zanardi , the essence of which 
we now summarize. Given some quantum system (com- 
puter) Q with the associated finite-dimensional Hilbert 
space 7i, we consider the error channel T with Kraus op- 
erators Ki . We define the interaction algebra & of T as a 
*-algebra generated by Ki. It is obvious that K, is an alge- 
bra with identity because of the condition J^i K*Ki = I- 
However, since the Kraus representation of a channel T is 
not unique, we must make sure that, for any two choices 
{Ki} and {K^} of Kraus representations of T, the corre- 
sponding interaction algebras are equal. Using the fact 
that any two Kraus representations of a channel are con- 
nected via 

Ki = J2 v ^ K ^ (65) 

where Vi^ are the entries of a matrix V with V*V = 1 
(assuming that one of the sets {Ki} and {K^} is padded 
with zero operators in order to ensure that they have the 
same cardinality), we see that it is indeed the case that 
the interaction algebra of a channel T does not depend 
on the particular choice of the Kraus operators. 

The existence of noiseless subsystems of Q with respect 
to T hinges on the reducibility of the interaction algebra 
M,. Since M. is a uniformly closed *-subalgebra of B(H), 
it is a finite-dimensional C*-algebra, and is therefore iso- 
morphic to a direct sum of r full matrix algebras, each of 
which appears with multiplicity rrii and has dimension nf 
(i.e., it is an algebra of rii x m complex matrices). Thus 
dim^ = Yli=i n l- The commutant &' of 8. is defined as 
the set of all operators X £ B(Tt) that commute with 
all K £ From the Wedderburn theorem pa] it follows 
that each K £ & has the form 

r 

K = l mi ® Ki , K t £ M ni (C) , (66) 
i=i 

and that each K' £ 8! has the form 

r 

K ' = K[ ® l ni , Ki £ M m% (C) . (67) 

i=l 

Thus dim&' = X)I=i m i- We have the corresponding 
isomorphism 

r 

H ~ C m ' ig) C ni , (68) 

8=1 

and each factor C mi is referred to as a noiseless sub- 
system because it is effectively decoupled from the error 



channel T. It is rather obvious that, in order to be of 
any use, a noiseless subsystem must be nontrivial, i.e., at 
least two-dimensional. Now, if the interaction algebra & 
is irreducible, then dim^' = 1, and no noiseless subsys- 
tems exist. There is a simple criterion of irreducibility 
of an algebra, the Schur's lemma [Q, which states that 
a *-algebra 21 is irreducible if and only if its commutant 
21' consists of complex multiples of the identity. We are 
now ready to state two main results of this subsection. 
Proposition 2. Let T be a strictly contractive unital 
channel. Then T admits no noiseless subsystems. 
Proof. Let us pick a Kraus representation {Ki} of T, 
and let & be the corresponding interaction algebra. We 
observe that if any X £ B(H) belongs to .ft', then X is a 
fixed point of T on B(7i). Indeed, 

X £ Si T{X) = KXKt =X^ KiK* . (69) 

i i 

Since T is unital, £\ K,K* = I, and thus T(X) = X. 
Now, if X £ then X* £ &' as well, which implies that 
X x := (X + X*)/2 and X 2 := (X - X*)/2i belong to 
&' . Therefore we only need to show that any selfadjoint 
X £ has the form rl for some r £ R. For any self- 
adjoint X, the operator \X\ :— (X 2 ) 1 / 2 belongs to the 
algebra generated by X 2 , so 

X = X* £# =^X ± := eJ*'. (70) 

Since X = X + — X_ and X± > 0, we reduce our task 
to proving that any positive X in 8! is a multiple of the 
identity. Without loss of generality, we may assume that 
\\X\\ X — 1. Since X > 0, we must have X £ S(Ti); since 
X belongs to the commutant of &, it is also a fixed point 
of T. Thus X = 1/ dim7i, and the commutant $! of the 
interaction algebra & consists of complex multiples of the 
identity. ■ 
The proof of Proposition 2 depends in an essential way 
on the uniqueness of the fixed point of a strictly contrac- 
tive channel, as well as on the condition satisfied by the 
Kraus operators of a unital channel. It turns out, how- 
ever, that the statement of Proposition 2 can be strength- 
ened to include all strictly contractive channels. 
Proposition 3. Let T be a strictly contractive channel. 
Then T admits no noiseless subsystems. 
Proof. Let A be the interaction algebra of the chan- 
nel T. Let us suppose, contrary to the statement of the 
Proposition, that T admits at least one noiseless subsys- 
tem (i.e., & is reducible). That is, there exists at least 
one j £ {1, . . . ,r} such that rrij,nj > 2 in Eqs. (|6^)-(|68|). 
Let fC be some closed subspace of TL. Restricting the 
channel T to the set 

S(K) := {p £ S(H) | supp p C JC} (71) 

(where supp p is the orthogonal complement of ker p), we 
note that, by definition, the contractivity modulus of the 
restricted channel cannot exceed the contractivity mod- 
ulus of T. Let Hj be the jth direct summand C m J <8> O 
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in Eq. 



Define the channel Ta as the restriction of T 



to S(TLj). Then any Kraus operator of Tj has the form 
Imj ® ^ where X M G M„ s (C) and 



E 



K*K 



I„ 



(72) 



Furthermore /c(Tj) < «(T) < 1. Now is the channel of 
the form id ®Sj, where Sj is the channel on S(C nj ) with 
Kraus operators Kn. As we have shown in subsection 



[II A, channels of this form are not strictly contractive 



(or even weakly contractive). Thus n{Tj) = 1, and the 
Proposition is proved, reductio ad absurdum. I 
The result of Proposition 3 is quite shocking as it 
unequivocally rules out the existence of noiseless sub- 
systems for any strictly contractive channel. From the 
standpoint of foundations of quantum theory, the im- 
portance of Proposition 3 lies in the fact that it estab- 
lishes nonexistence of noiseless subsystems for a wide 
class of physically realizable quantum computers on the 
basis of a minimal set of assumptions. Furthermore, from 
the mathematical point of view, it is rather remarkable 
that strict contractivity of a channel already implies irre- 
ducibility of its interaction algebra. We must, however, 
hasten to emphasize that, despite its sweeping general- 
ity, Proposition 3 should not be considered as a proof of 
impossibility of building a reliable quantum computer. It 
merely rules out the possibility of building quantum com- 
puters with perfect protection against errors modelled by 
strictly contractive channels. 



D. Approximate error correction 

At this point we must realize that the results of the pre- 
vious subsection are not as unexpected as they may seem. 
After all, nothing is perfect in the real world! Therefore, 
our error-correction schemes must, at best, come as close 
as possible to the perfect scenario. Of course, the precise 
criteria for determining how close a given error-correction 
scheme is to the "perfect case" will vary depending on the 
particular situation, but we can state perhaps the most 
obvious criterion in terms of distinguishability of chan- 
nels. 

Let us first phrase everything in abstract terms. Let 
the error mechanism affecting the computer be modelled 
by some channel T. We assume that there exists some 
positive 8 < 1 which, in some way, characterizes the 
channel T (it could be given, e.g., by the minimum of 
the operator norms of the Kraus operators of T, and 
thus quantify the "smallest" probability of an error oc- 
curring). Let H be the Hilbert space associated with 
the computer. Then, for each e > 0, we define a (e, 8)- 
approximate error- correcting scheme for T to consist of 
the following objects: 

(1) an integer n > 1, 

(2) a Hilbert space 7i e xt with dim7i oxt > dim7i, 



(3) a channel E : S(H) -> 5(W ext ), 

(4) a channel T : S(7i xt) — > S(H ex t), and 

(5) a completely positive map T corl : iS(TYcxt) * 

such that the channel T depends uniquely on n, Ti ex t, T, 
and E; the CP map T corr is correctable (say, in the Knill- 
Laflamme sense, or through other means, depending on 
the particular situation); and we have the estimate 



T - T„ 



< S n < e. 



cb 



(73) 



Let us give a concrete example in order to illustrate 
the above definition. Suppose that the channel T is of 
the form id+S* with ||5|| cb < 5. Then, for any n, we can 
write 



T 



id- 



E 

AC{1 «} k=l 

0<|A|<7i 



(74) 



where \A\ denotes the cardinality of the set A, and 
la '■ {1, • • • , n} — > {0, 1} is the indicator function of A. 
We use the convention that, for any map M, M° = id. 
In other words, the summation on the right-hand side 
of Eq. (|7^ ) consists of tensor product terms with one 
or more identity factors. For the last term, we have 
\\S® n \\ ch <S n . 

In this case, given some e > 0, we pick such n that 
8 n < e and let H cxt := H® n . If the CP map given by 
the sum of the first two terms on the right-hand side of 
Eq.(f74j) is correctable on some subspace K, of H ex t, then 
the channel E is defined in a natural way through the 
composition of the following two operations: (a) adjoin- 
ing additional n — 1 copies of Ti, each in some suitable 
state po-, and (b) restricting to the subspace /C. This way, 
we obviously have T := T®" and 



E 



(75) 



AC{l,...,n} fc= 
0<|A|<n 



The estimate (|§ holds because T - T corr = S® n . We 
note that this construction results in a quantum error- 
correcting code that corrects any n — 1 errors. We can 
use similar reasoning to describe quantum codes that cor- 
rect k < n errors. 

Constructing 7i cx t as a tensor product of a number 
of copies of Ti, the Hilbert space of the computer, evi- 
dently leads to the usual schemes for fault-tolerant quan- 
tum computation |36|. Other solutions, such as embed- 
ding the finite-dimensional Hilbert space TL in a suitable 
infinite-dimensional Hilbert space (e.g., encoding a qubit 
in a harmonic oscillator |37|]), can also be formulated in 
a manner consistent with our definition above. 
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Let us now address approximate co rrecta bility of 
strictly contractive errors. In subsection [V E we have 



demonstrated that, in the absence of error correction, the 
sensitivity of quantum memories and computers to such 
errors grows exponentially with storage and computation 
time respectively. Let T be a strictly contractive error 
channel. It is obvious that the appropriate approximate 
error-correction scheme must be such that the contrac- 
tion rate of the "encoded" computer, where the errors 
are now modelled by the channel T, is effectively slowed 
down. In some cases, straightforward tensor-product re- 
alization may prove useful (e.g., when the product chan- 
nel T <S> T is not strictly contractive) . We must recall 
that, for any channel S, a necessary condition for cor- 
rectability is k(S) = 1. Thus, if we can find a suitable ap- 
proximate error-correcting scheme where T would be well 
approximated by some channel T corr with k(T coii ) = 1, 
we may effectively slow down the contraction rate by pro- 
tecting the encoded computer against errors modelled 
by Tcorr- A more ingenious approach may call for re- 
placing circuit-based quantum computation with that in 
massively parallel arrays of interacting particles; several 
such implementations have already been proposed [j38| . 
It is quite likely that the possible "encodings" of quan- 
tum computation in these massively parallel systems [|3Q| ] 
may offer a more efficient implementation of approximate 
error correction. 

Finally, we should mention that the idea of "approx- 
imate" noiseless subsystems has already been explored 
by Bacon et al. [Q. In their work, it is argued that 
the symmetry, which is required of a channel in order 
for noiseless subsystems to exist, is generally broken by 
perturbing the channel. They show that, if the perturba- 
tions of the channel are "reasonable," then the noiseless 
subsytem is stable to second order in time. We must 
reiterate that the negative results we have stated in the 
previous subsections refer only to nonexistence of "per- 
fectly" noiseless subsystems; in the real world, we would 
have no choice but to settle for " almost perfect" anyway. 



V. CONCLUSION 

In this paper, we have offered an argument that er- 
rors in physically realizable quantum computers are natu- 
rally modelled by strictly contractive channels, i.e., chan- 
nels that uniformly shrink, in the trace norm, the set of 
all density operators of the system under consideration. 
In particular, no two density operators in the image of 
a strictly contractive channel have orthogonal supports, 
which implies that any measurement designed to distin- 
guish between these density operators will err with prob- 
ability bounded away from zero. This implies, in turn, 
that there exists some precision threshold e > such that 
any two density operators p, a with ||p — °"Hi < e cannot 
be distinguished by a particular experimentally available 
measuring apparatus. 



We can turn this reasoning around by first postulating 
the existence of a precision threshold e that would quan- 
tify resolving power of the least precise instrument em- 
ployed in the experiment. As we have argued, the phys- 
ical interpretation of the precision threshold boils down 
to limits on our ability to distinguish between density 
operators. A nonzero lower bound on the probability of 
error in optimum quantum hypothesis testing can thus be 
taken as an indication that the combined influence of en- 
vironmental noise and experimental imprecisions (which, 
in fact, are quite likely to be caused by indelible quantum- 
mechanical effects, such as vacuum fluctuations) can be 
economically captured by the concept of a strictly con- 
tractive channel. 

As we have shown, the set C sc of all strictly contrac- 
tive channels on a given system Q is dense in the set C 
of all channels on Q. Since no finite-precision measure- 
ment will be able to distinguish between an arbitrary 
channel T and some strictlly contractive channel T', it 
is reasonable to ascribe to strictly contractive channels 
the property of "experimental reality," just as we would 
ascribe this property to elements of the set [0, 2ir) n Q 
(where Q is the set of rational numbers) in any experi- 
ment involving finite- precision measurements of angles. 

In light of this interpretation, it is important to in- 
vestigate the robustness of quantum memories and com- 
puters in the presence of strictly contractive errors. We 
have found that, in the absence of error correction, any 
state stored in a noisy quantum register converges expo- 
nentially fast to the fixed point of the error channel T, 
and the rate of convergence is independent of dimension 
of the register Hilbert space. In other words, sensitivity 
of quantum registers to strictly contractive noise is an 
intensive property, i.e., independent of the register's size. 
Similarly, computations performed on a noisy quantum 
computer with different initial states quickly yield indis- 
tinguishable results, again at a rate that does not depend 
on the computer's size (number of qubits). Furthermore, 
the property of strict contractivity turns out to be strong 
enough to proscribe the existence of noiseless subsystems 
of the computer affected by any strictly contractive error 
channel. 

However, these results are more of a blessing than 
a curse for the future of quantum information process- 
ing: they certainly indicate that the successful solution 
of problems faced by researchers in this field will require 
models of computers far more ingenious than networks of 
one- and two-qubit gates. As we have mentioned above, 
systems of interacting particles (or quantum cellular au- 
tomata) may well prove to be a viable medium for exper- 
imental realization of large-scale quantum computers. In 
this respect, we would like to point out a possible con- 
nection between strictly contractive channels and ergodic 
quantum cellular automata [fhl| . A cellular automaton is 
ergodic if it possesses a unique invariant state which it 
reaches irrespective of initial conditions, and this is ex- 
actly the property shared by quantum systems under the 
influence of strictly contractive errors. As an example, 
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let us consider information storage in a quantum cellu- 
lar automaton. It is essential that this automaton be 
nonergodic, for otherwise it would not be able to "re- 
member" anything. Assuming that each cell (site) of the 
automaton is under the influence of some strictly con- 
tractive error channel T, an interesting problem would 
be to devise such a transition rule that the automaton 
would not be ergodic. In this respect, we should mention 
that, while T is a strictly contractive channel, it is not at 
all obvious whether T ® T is strictly contractive as well: 
it has a unique fixed point among the product density 
operators, namely pr ® Pt, but there may also be an- 
other fixed point of T ® T that is not a product density 
operator. 

Finally we mention one more point worth exploring. In 
our discussion of physically realizable quantum comput- 
ers, we have implicitly assumed that the (im)precision 
of all experimentally available procedures can be traced 
back to the (im)precision of state preparation, quanti- 
fied by some threshold value e (i.e., when we say that 
state p has been prepared, we mean that any state a 
with \\p — <j\\ x < e may have emerged from our preparing 
apparatus), as well as the (im) precision of measurements 
(we would not be able to distinguish any two states p, a 
with Hp — cr 1 1 x < e). This is suggestive of Ludwig's ax- 
iomatics of quantum theory ]4^ ] , and it would be theoreti- 
cally rewarding to consider physically realizable quantum 
computers from this axiomatic perspective as well. 



ACKNOWLEDGMENTS 



The author would like to thank: G.M. D'Ariano, M. 
Ozawa, and H.P. Yuen for enlightening discussions and 
constructive criticism; D.A. Lidar for interesting com- 
ments and pointers to useful references; and the ref- 
eree for helpful suggestions. This work was supported 
by the U.S. Army Research Office through MURI grant 
DAAD19-00-1-0177. 



Electronic address: maxim@northwestern.edu; phone: 

(847) 491-3063; fax: (847) 491-4455. 
[1] W.G. Unruh, Phys. Rev. A 51, 992 (1995); G.M. Palma, 

K.-A. Suominen, and A.K. Ekert, Proc. R. Soc. London, 

Ser. A 452, 567 (1996). 
[2] From this point on, we will refer to quantum computers 

as simply "computers" for brevity. 
[3] E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 
[4] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 

(1997). 

[5] D.A. Lidar, IL. Chuang, and K.B. Whaley, Phys. Rev. 
Lett. 81, 2594 (1998). 



[6] E. Knill, R. Laflamme, and L. Viola, Phys. Rev. Lett. 
84, 2525 (2000). 

[7] P. Zanardi, Phys. Rev. A 63, 012301 (2001). 

[8] Please note that we must be careful to distinguish preci- 
sion, which is a figure of merit for individual operations, 
from accuracy, which is a term that refers to repeated 
operations of the same kind (cf. @ for a thorough dis- 
cussion). In our case, precision of a physically realizable 
computer is a measure of the closeness of its output state 
to that of the ideal computer, while accuracy quantifies 
the closeness of output states, generated by the computer 
on repeated runs, in terms of the closeness of the corre- 
sponding input states. 

[9] T. Brody, The Philosophy behind Physics (Springer- 
Verlag, Berlin, 1993), pp. 139-156. 
[10] R. Durrett, Probability: Theory and Examples, 2nd ed. 

(Wadsworth, Belmont, 1996). 
[11] E. Bernstein and U. Vazirani, SIAM J. Comput. 26, 1411 
(1997). 

[12] K. Kraus, States, Effects, and Operations (Springer- 
Verlag, Berlin, 1983). 

[13] D.A. Meyer, Phys. Rev. Lett. 83, 3751 (1999). 

[14] O. Bratteli and D.W. Robinson, Operator Algebras and 
Quantum Statistical Mechanics, vol. 1, 2nd ed. (Springer- 
Verlag, 1987), pp. 42-64. 

[15] D. P. Zhelobenko, Compact Lie Groups and Their Rep- 
resentations (Amer. Math. Soc, Providence, 1978). 

[16] R. Bhatia, Matrix Analysis ( Springer- Verlag, New York, 
1997), ch. IV. 

[17] A. Yu. Kitaev, Russian Math. Surveys 52, 1191 (1997). 
[18] A.S. Holevo, Probabilistic and Statistical Aspects of 

Quantum Theory (North-Holland, Amsterdam, 1982), p. 

28. 

[19] V.I. Paulsen, Completely Bounded Maps and Dilations 

(Longman Scientific & Technical, Harlow UK, 1986). 
[20] D. Aharonov, AYu. Kitaev, and N. Nisan, in Proc. 30th 

ACM Symposium on Theory of Computing, p. 20, 1997. 
[21] G. Giedke, H.J. Briegel, J.I. Cirac, and P. Zoller, Phys. 

Rev. A 59, 2641 (1999). 
[22] A.S. Holevo, R.F. Werner Phys. Rev. A 63, 032312 

(2001). 

[23] M. Raginsky, "A fidelity measure for quantum chan- 
ne ls" (Phys. Let t. A, in press); LANL preprint 



[24 
[25 

[26 
[27 
[28 
[29 
[30 
[31 



quant -ph/010710£ . 

G.M. D'ariano, P. Lo Presti, and M.G.A. Paris, LANL 



preprint |quant-ph/0109040 . 

M. Reed and B. Simon, Methods of Modern Mathemati- 
cal Physics I: Functional Analysis (Academic Press, San 
Diego, 1980), ch. V. 

I. Stakgold, Green's Functions and Boundary Value Prob- 
lems, 2nd ed. (Wiley, New York, 1998), ch. 4. 
G. Mahler and V. Weberruss, Quantum Networks, 2nd 
ed. (Springer- Verlag, 1998), ch. 2. 

D. Bruss, L. Faoro, C. Macchiavello, and G.M. Palma, J. 
Mod. Opt. 47, 325 (2000). 

M.B. Ruskai, S. Szarek, and E. Werner, LANL preprint 



quant -ph/0101003 



GG Amosov, A.S. Holevo, and R.F. Werner, Probl. In- 
form. Transm., 36, 25 (2000). 

C.W. Helstrom, Quantum Detection and Estimation 



13 



Theory (Academic Press, N ew York, 1976), ch. IV. 
[32] H.P. Yuen, LANL preprint |quant-ph/0006109| v7 (Ap- 
pendix A). 

[33] Our definition of a quantum circuit differs from the con- [39] 
ventional one. See, e.g., A.C.-C. Yao, in Proc. 34th IEEE 
Symposium on Foundations of Comp. Sci., p. 352, 1993. 

[34] A. Barenco, CH. Bennett, R. Cleve, D.P. DiVincenzo, 
N. Margolus, P. Shor, T. Sleator, J.A. Smolin, and H. 
Weinfurter Phys. Rev. A 52, (1995). [40] 

[35] L. K. Grover, Phys. Rev. Lett. 79, 325 (1997). 

[36] J. Preskill, Proc. R. Soc. London, Ser. A 454, 385 (1998). [41] 

[37] D. Gottesman, A.Yu. Kitaev, and J. Preskill, Phys. Rev. [42] 
A 64, 012310 (2001). 

[38] P. Zanardi and F. Rossi, Phys. Rev. B 59, 8170 (1999); 



H.J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 
910 (2001); R. Raussendorf and H.J. Briegel, Phys. Rev. 
Lett. 86, 5188 (2001). 

The term "massively parallel" here refers to the symmet- 
rical nature of quantum arrays of interacting particles. 
This usage is different from the usual connotation of mas- 
sive parallelism, intrinsic in the mathematical structure 
of quantum theory. 

D. Bacon, D.A. Lidar, and K.B. Whaley, Phys. Rev. A 
60, 1944 (1999). 

S. Richter and R.F. Werner, J. Stat. Phys. 82, 963 (1996). 
G. Ludwig, An Axiomatic Basis for Quantum Mechanics, 
vols. 1 and 2 (Springer- Verlag, Berlin, 1985). 



14 



